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Common beliefs about unitarity are not reliable, and we do not know how to 
apply DGLAP evolution at small x. Together with the big discrepancy between 
the measurements of the total cross section at the Tevatron, a consequence is that 
the cross section at the LHC could be anywhere between 90 and 160 mb. 

PACS numbers: 11.55.Jy 13.60.Hb 13.85.Dz 13.85.Lg 

1. Introduction 

The LHC will probe physics under extreme conditions. This means that past 
physical intuition may be unreliable, and models that depend on intuition 
may break down. Therefore it is important to ask which of our current 
beliefs have a sound theoretical basis, and which are just based on folklore. 

I will concentrate on two topics where our fundamental understanding is 
particularly uncertain: 

• unitarity 

• DGLAP evolution at small x. 

I will show that our lack of understanding of these has serious consequences 
for what is probably the first thing that will be measured at the LHC, the 
total cross section. The best estimate has a huge error: 

,LHC 



^LHc ^ ± 35 mb (1.1) 



My understanding of most of the material that I review here derives from 
my work over the years with Sandy Donnachie. Further details may be 
found in our book[lJ. 

* Presented at Epiphany meeting to commemorate Jan Kwiecinski, Krakow, January 
2009 
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2. Unitarity 

For an elastic hadron-scattering amplitude, the unitarity equation reads 
^C_J— ^ — C*)— ^ ' — inelastic terms 

or 

Im a^(s) = |a^(s)p + inelastic terms (2-1) 
so that the partial-wave amplitude obeys 

\ae{s)\ < 1 (2.2) 
A well-known consequence of (12. 2p is the Froissart-Lukaszuk-Martin bound [2] [3]: 



a™^(.) < ^ log\s/so) (2.3) 



7r 



At LHC energies, for reasonable values of the unkown scale sq, this gives a 
bound of several barns, and so it is not a useful constraint. 

A more useful bound is that of Pumplin[3]: 

^ELASTIC ^ 1 ^TOTAL (2.4) 

The exchange of a single pomeron exchange JP gives 

^TOTAL ^^e ^ ^ Q Qg 

and 

ELASTIC 

s"^ (2.5) 



ELASTIC 
OCT 2e 
r 

dt 



which clearly causes this bound to be violated at large enough s. The 
remedy is to sum single- iP, double- iP, . . . exchanges: 
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However, even though we started trying to learn how to do this more than 
40 years ago, we still do not know how. 

The best we can do, although it is certainly wrong[lJ, is to use an eikonal 
formalism. Write the amplitude as a 2-dimensional Fourier integral 

A{s, -q2) =Aj (fhe-''^-^A{s, b^) (2.6) 

Define ^(•s, b) = — log(l + 2iA/s) so that 

i(s,b2) = iis(l-e-^('''')) (2.7) 

Then the bound (j2.2p on the partial-wave amplitude may be shown to be 
equivalent to 

Rex(s,6)>0 (2.8) 

which is easy to impose. If we expand the exponential in (|2.7|) in powers of 
X we have 

A(s, -q2) = 2is J d^e-'^^-^il - g-^^"'^)) 

= 2i.|d26e-^-^(x-|- + |-...) (2.9) 

So far, the equations are certainly correct. But we do not know what to take 
for xis,b). An obvious choice is to approximate it by single-iP exchange. 
Although we do not know how to calculate double- iP exchange, we do know 
something about its general structure, and the second term in (j2.9|) has 
the right structure. Similarly, the third term has the right structure to 
represent triple-iP exchange, and so on. Nevertheless, this is wrong, for 
various reasons [1]. One one is that double- iP exchange obviously depends 
on two-quark correlations in the proton wave function, which are absent in 
this procedure. 

There has been a lot of talk about the consequences of unitarity for processes 
with leptons or photons in the initial state; for example it is believed to 
lead to what is known as saturation. It is important to understand that 
basic theory alone does not allow one to conclude anything, without feeding 
in extra assumptions. That is, there may well be no unitarity bound on 
F2{x,Q'^). The reason is that the analogue of (|2.ip for j*p scattering 




2 

+ inelastic terms 



4 



KRAKOW'PUB PRINTED ON MARCH 9, 2009 




1e-06 1e-05 0.0001 0.001 



X 

Fig. 1. Simple fit (|3.4p to data for F2{x, Q^) at small x 



is not true. This is the case even for ^yp scattering: for ah we know, the 
■jp total cross section could continue to increase indefinitely with increasing 
energy. 



3. DGLAP evolution at small x 

Claims by experimentalists and theorists to be able to extract parton den- 
sities with high accuracy from data, for use at the LHC, must be regarded 
with caution. This is because we do not know how to handle DGLAP evo- 
lution at small x. 
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The singlet DGLAP eqation is: 



(3.1) 

The terms of the perturbation expansion of 'P{z, as{Q'^)) diverge like 1/ z at 
z = 0. For small x the integration extends to small values of z and therefore 
it is wrong to expand 'P{z, as{Q'^)) in powers of as- 

A possible exception is when u(x, Q^) rises steeply with To understand 
this, approximate u(x, Q^) in some range of x at some value of by 

and insert this into the DGLAP equation (j3.ip . This gives 

d_ 

di 



logf(Q2) = p(iV = e,a,(Q2))_ f dzz'V{z,as{Q^)) (3.3) 

JO 



with P the Mellin transform of P. The last term ~ x'^ and so is negligible 
at small x if e is some way above 0. A pole of P(z, agiQ"^)) at z = reflects 
itself in a pole of P(A^, as((5^) at = 0. So if e is some way above we are 
at some distance from this pole and then also it should be safe to expand 
the first term in powers of a^. By doing so, and dropping the last term in 
(j3.3p . we obtain a simple differential equation for f((5^)- 

The simplest fit to F2 at small x is a combination of two powers of x, 
hard-pomeron and soft-pomeron: 

F2{x. Q') = /o(Q')x-^° + /i(Q')x-^i 

eo ~ 0.4 ei = 0.0808 (3.4) 
The fit has just 5 free parameters, including eo. See figure [TJ 

Donnachie and I made this fit[5j purely phenomenogically, but we then 
found[7] that its output for the coefficient function fo{Q'^) in (j3.4p obeys 
the DGLAP evolution differential equation to very high accuracy, both at 
LO and at NLO. See figure [2j DGLAP evolution is supposed to be valid 
only for large and we can see from the figure that this means at least 5 
to 10 GeV^. It certainly does not make sense to use it down to 1 or 2 GeV^, 
as is often done. 

As I have explained, DGLAP evolution breaks down for the soft-pomeron 
coefficient function fi{Q'^), as ei is too close to 0. The conventional approach 
to DGLAP evolution, used by many theorists and experimentalists to ex- 
tract what are claimed to be highly accurate parton distributions, amounts 
to ignoring this difficulty. 
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Our simple fit (j3.4|) applies only at small x, which is why figure [T] shows 
data only for x < 10^^. If we want to extend it to larger values of x, we 
should add in a Regge term x"*^^ with ^ — |, corresponding to /2 and 
02 exchange. Also, the simple powers of x in (j3.4p must be multiplied by 
functions that go to as x ^ 1. We do not know what these should be, so 
Donnachie and I took[5j just the powers of (1 — x) given by the dimensional 
counting rules. This is certainly too simple to be correct but is better than 
doing nothing. It is also astonishingly successful: see figure [3l which shows 
also what the fit gives for the ^yp total cross section. Compared with figure 
13.41 only 2 free parameters have been added, both for the Regge term. 



4. Total cross section at the LHC 

Given that the data for F2{x,Q'^) respond so well to a fit that includes a 
hard pomeron, it is natural [6j to include such a term also in fits to pp and pp 
scattering. That is, for each of o-{pp), a{pp), a{'yp) include hard pomeron, 
soft pomeron and Regge exchange: 

a = Xos« + Xis'' + Xrs'^ (4.1) 

The result is shown in figured! In the left-hand figure, the lowest curve is the 
hard-pomeron contribution. The right-hand figure shows the extrapolation 
to LHC energy of the two fits, this new one and the old fit without a hard 
pomeron term. 

Notice the familiar and long-standing discrepancy between the E710 and 
CDF measurements at the Tevatron. If the upper CDF measurement should 
be correct, it surely is a sign that something new is beginning to become 
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Fig. 3. The fit to F2 extended to larger values of x, and to real-photon data 
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Fig. 4. Fits to pp and pp total cross sections 



important at Tevatron energy, with the consequence that the LHC cross 
section wih be large. This observation does not depend on any particular 
theoretical explanation of the data. 

However, the large prediction for the LHC total cross section when the hard 
pomeron is included will lead some to worry about unitarity. Some even 
worry about unitarity for the lower curve. Donnachie and I stressed when we 
made the original DL fit [9] that the power ei ~ 0.08 was an effective power 
that already, to some extent, includes unitarity corrections; nevertheless 
Alan Martin and collaborators believe [8] that this is not enough to take 
account of unitarity and so predict that the LHC cross section will be about 
90 mb. 

I have explained that nobody knows how to calculate the effects of unitarity. 
I will now describe an attempt to do so, which should not be taken at all 
seriously. 

Consider pp and pp elastic scattering. At existing energies, soft pomeron 
exchange dominates and its contribution has been known for nearly 40 years: 

da [3/3iFi(t)]4 



dt 4tt 



-(a'is)2(^i+"i*) (4.2) 



Here Fi{t) is the elastic form factor of the proton and /3i and ei are known 
from a'^'^'^. As long ago as 1973, my then student Jaroskiewicz[10j found 
that data fix the only free parameter a'^ to be 0.25 GeV~^. This is done 
by using very-small-t data at some energy: see the upper part of figure ([5]). 
Then the formula (j4.2p fits the data out to rather larger values of t at that 
energy, as is seen in in the lower plot in the figure. Because Fi (t) is raised to 
the 4th power the fit is rather sensitive to it. I do not understand why the 
proton's electromagnetic form factor should be appropriate, since pomeron 
exchange has the opposite C parity from photon exchange. 
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Fig. 5. pp elastic scattering data at = 53 GeV. The curves do not include 
photon exchange 



As is well known, at small t pp elastic scattering displays shrinkage: da/dt 
becomes steeper as the energy increases. The rate of shrinkage is determined 
by the value of a'l, and 0.25 GeV~^ describes the data very well at all the 
different energies that have been measured [T]. 

At large values of greater than about 3 GeV^, the data take on a different 
character. They are described very well by 



da 

— = 0.09 r 

dt 



(4.3) 



and are independent of energy. See figure [H The form (j4.3p is what is 
obtained from the triple-gluon-exchange mechanism of figure [71 This raises 
an interesting question pT]: what if one replaces each gluon with a hard 
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Fig. 6. Large-i data for pp elastic scattering at various energies, with the fit (|4.3p 




Fig. 7. Triple-gluon exchange 



pomeron? This might provide a mechanism which, while it is too small to 
be seen at existing energies, grows rapidly with energy and is large at LHC 
energy. That is, da/dt at the LHC might be large at large t. 

At energies corresponding to the data plotted in figures [5] and [6l pp elastic 
scattering data display a striking dip structure at values of t in between 
those of the two figures. See figure El It is not easy to generate a dip: 
the real and imaginary parts of the amplitude must be very small at the 
same t, which requires something of an accident. It needs at least three 
contributions: probably iP, IPIP and ggg. However, the third of these 
three terms has the opposite C parity from the first two, which led us to 
predict [12] that pp scattering should not have a dip. As the figure shows, 
this was later confirmed. An exchange such as ggg which has negative C 
parity is known as odderon exchange; it is a mystery why odderon exchange 
has not been detected at smaller values of t. 
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Fig. 8. pp and pp elastic scattering data. The left hand figure is for pp and the 62 
GeV data are multiplied by 10. The data in the right-hand figure is at 53 GeV and 
the upper points are pp, the lower pp. 




Fig. 9. pp elastic scattering data from the Tevatron; the curve corresponds to the 
exchange of single hard and soft pomcrons 



The fit shown in figure[5]to elastic scattering at relatively low energy remains 
good at very small t as the energy is increased but becomes less good at 
larger values of t. Figure [9] shows data for pp elastic scattering from the two 
Tevatron experiments. The curve corresonds to the exchange of a single soft 
pomeron plus that of a single hard pomeron. I have explained that we do not 
know how to calculate the exchange of two pomerons, but we know enough 
about its general features to know that including this will pull da/dt down 
at larger t. As a very crude model that includes the two-pomeron exchanges 
FqIPo, IPiJPi and FqPi (where again the subscripts and 1 denote the hard 
and soft pomerons) , I calculated the 6-space amplitude for the sum of the JPq 
and iPi exchanges, and squared it to simulate the sum of the IP IP exchanges: 



(4.4) 
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Fig. 10. pp and pp elastic scattering data at various with fit JP + PIP + ggg 
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Fig. 11. Extrapolations to LHC energy of fits to the total cross section. The upper 
curve corresponds to the fits to the amplitude shown in figure [TUl including PIP 
exchanges, while the lower curve omits any hard-pomeron contribution 



I chose the value of A to cancel the imaginary part of the amplitude so as 
to get the dips in pp scattering at the right value of t, and added in ggg 
exchange. The large-t form ([6]) of the latter needs to be modified at smaller 
t so that it does not diverge, and one does not know how to do this, but I 
guessed a form with a single parameter which I chose so as the optimise the 
fit to the dip structure. My best fit to the data is shown in figure [TOl 

The result for the total cross section at the LHC is that including the two- 
pomeron exchanges pulls the prediction down from 160 to 125 mb: see 



krakow'pub printed on March 9, 2009 



13 



figure [TH The lower curve in the figure corresponds to single soft pomeron 
exchange only. So, since my attempt to include the double exchanges is 
surely very crude, I have to conclude that the LHC total cross section could 
be anywhere between 100 and 160 mb. Remember, though, that it might 
be even smaller [8]. 

5. Summary 

• CT^"^ = 125± 35 mb 

• We do not know how usefully to impose unitarity — eikonal-type models 
are surely too simple 

• Unitarity does not constrain lepton or photon-induced cross sections 

• We still cannot calculate JPIP exchange - even after more than 45 years 

• There are severe mathematical problems with DGLAP at small x 

• DGLAP cannot be used below = 5 GeV^ 

• Regge fits to F2{x^ Q^) are the simplest — and probably the most correct 

• Elastic scattering at the LHC at large t may be surprisingly large 

• If the CDF Tevatron cross section is correct, something dramatic must 
happen — independently of any theory! 
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